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1. Introduction 



Coherent neuronal oscillations have been implicated widely as a correlate of brain function (Gray, 
1994; Llinas and Ribary, 1993). However, our understanding of the mechanisms underlying such 
activity is still in its infancy. (See Jefferys et al. (1996) and Ritz and Sejnowski (1997) for recent 
reviews). It has been demonstrated experimentally (Whittington et al, 1995; Jefferys et a/., 1996), 
computationally (Wang and Rinzel, 1992; Whittington et al., 1995), and analytically (van Vreeswijk 
et al., 1994; Gerstner et al., 1996; Terman et al., 1997) that mutual inhibition can generate stable 
synchronous oscillations. Previous analytical work mostly concentrated on the mechanisms respon- 
sible for synchronization. Here we address the mechanisms for controlling the network frequency 
and the interplay between frequency and synchronization. 

Previously (White et al., 1998), we showed that synchronization of inhibitory networks can 
occur in a wide range of frequencies for homogeneous neurons, but in the presence of heterogeneity 
network synchronization is very fragile in some parameter regimes. We have identified a pair of 
parameter regimes, denoted 'phasic' and 'tonic'; in the former some coherence is maintained in the 
presence of mild heterogeneity, while in the latter it is lost. In this paper, we relate these regimes 
to the parameters that determine the frequency of the network when it is synchronized. 

There are three important time constants in the network dynamics. One is an intrinsic time con- 
stant of the uncoupled cell, dependent on the conductances and the capacitance of the membrane. 
The second is the decay time of the inhibition. The third, the network period, depends on the other 
two time scales and other parameters, notably the synaptic conductance and injected current. The 
aim of this paper is to understand how the first and second time constants influence the third. 

We show that there are three asymptotic regimes, in which the network period varies differently 
as the inhibitory time constant is changed. Two of these are the tonic and phasic (White et al., 
1998) regimes mentioned above and the third we call 'fast'. In the tonic regime, the period is 
small compared with both the intrinsic time scale and the synaptic decay time. In this regime, we 
show that the network frequency is only weakly dependent on the inhibitory decay time; indeed, 
the frequency is affected mainly by the average amount of inhibition, not by the time course of 
that inhibition. In the phasic regime, the intrinsic membrane time scale is small compared with 
the network period and the synaptic decay time. In this regime, we show that the network period 
is proportional to the decay time, with the proportionality constant a function of other network 
parameters. In the fast regime, the synaptic decay time is short compared to the network period. 
In this regime, the period is dominated by the intrinsic time constant. We give constraints on the 
network parameters (intrinsic membrane time scale, synaptic decay time, synaptic conductance and 
applied current) for the network to fall into each of the three asymptotic regimes. 

The three regimes can be related to the behavior of the network in the presence of mild het- 
erogeneity , as shown numerically in White et al. (1998). In the fast regime, the synaptic influence 
acts as a brief inhibitory pulse which has been shown to lead to stable anti-synchronous oscillations 
(Friesen, 1994; Perkel and Mulloney, 1974; Skinner et al., 1994; van Vreeswijk et al., 1994; Wang 
and Rinzel, 1992). In the phasic and tonic regimes, synchrony through inhibition is fairly robust for 
homogeneous networks of neurons over a wide range of network frequencies (van Vreeswijk et al., 
1994; Gerstner et al., 1996). However, it is very fragile when even mild heterogeneity is included 
(Golomb and Rinzel, 1993; Wang and Buzsaki, 1996; White et al., 1998). We showed that the loss of 
coherence of the network happens by different mechanisms in the tonic and phasic regimes (White 
et al., 1998). In the tonic regime, mildly heterogeneous networks are effectively de-coupled and 
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exhibit asynchrony (loss of phase relationships among the cells); in the phasic regime, the network 
can lose some coherence via suppression, in which the slower cells receive enough inhibition to 
prevent their firing. 

Our analytical work uses a reduced membrane description for the neuron with a biophysically- 
based synaptic model. The analysis is for a homogeneous network; when considering the frequency 
of a synchronous solution, we can regard the network as a single self-inhibited neuron. We show that 
in the phasic regime, the synaptic current dominates and the actual intrinsic membrane dynamics 
are not as important for the firing frequency. In the tonic regime, the intrinsic dynamics do become 
more important. We give a comparison between our analytical estimates for the period and those 
obtained numerically for conductance-based neuron models. 

The work in this paper, in conjunction with White et al. (1998) shows that, in the presence of 
mild heterogeneity, the control of frequency is strongly tied to the time scale of the mechanism 
that produces the synchronization. That is, in order to have coherence for a mildly heterogeneous 
system, the latter must be in the phasic regime, which then implies that the frequency is controlled 
by the time constant of the inhibitory decay. It was found experimentally and numerically that 
the frequency of coherently firing interneurons in the CAl region of the hippocampus is strongly 
dependent on the decay time of the inhibitory synapse (Whittington et al., 1995; Jefferys et al., 1996; 
Traub et al., 1996a; Wang and Buzsaki, 1996). The analytical work given in this paper clarifies the 
reasons for this frequency behavior as well as the remark found in many papers that the inhibitory 
decay time can be a critical factor in the determination of the network frequency (Destexhe et al., 
1993; Skinner et al, 1993; Kopeh and LeMasson, 1994). 



We consider the frequency control of so called Type I neuronal dynamics (Hansel et al., 1995; 
Ermentrout, 1996) and examine simplifications that would make it more amenable to analysis. These 
neurons are distinguished by the fact that they have positive phase response curves (PRC) (Hansel 
et al., 1995). i.e. a positive depolarizing current always advances the time of the next spike. It 
has been shown recently that neurons which admit very low frequency oscillations near the critical 
applied current are Type I (Ermentrout, 1996) although the converse is not necessarily true. We 
note that the type of bifurcation to firing of the neuron is not important for our analysis. Type 
I neurons have been used to represent inhibitory interneurons in the hippocampus (Traub et al., 
1996a; Wang and Buzsaki, 1996; White et al, 1998). 

We consider an inhibitory network of single-compartment neurons with membrane potential 
dynamics of the form (Rinzel and Ermentrout, 1989; Ermentrout, 1996; Hansel et al., 1995; Wang 
and Buzsaki, 1996; White et al, 1998) 



where /j is an applied current, /ion are the ionic currents responsible for spike generation and recov- 
ery and Is is the synaptic current induced by the spikes of other neurons coupled pre-synaptically. 
The ionic currents lion are functions of potential and the dependent dynamical variables, which are 
in turn governed by a system of differential equations. In the Appendix we give the equations for 
two conductance-based neuron models (White et al, 1998; Ermentrout and Kopell, 1998), which 
are Type I or at least very close to Type I (data not shown). Figure Q shows example voltage traces 



2. Neuron model 




(1) 
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for the two models. The couphng between the neurons is exclusively through chemical synapses 
represented by the synaptic current which takes the form Ig = gS{t){V — Vs), where 

^ = aF{Vpre){l-S)-pS, (2) 

a and /3 are respectively the synaptic rise and decay rates, Vpre is the membrane potential of the 
pre-synaptic neuron, and F{V) = 1/(1 + exp[— y]). We will often consider the synaptic decay time 

We approximate the full conductance-based dynamics with integrate-and-fire dynamics where 
the firing of the neuron is represented by the resetting of the membrane potential whenever it 
crosses a threshold. Our justification for using this approximation hinges on three considerations: 
1) both the conductance-based and integrate-and-fire models are Type I (in the sense of positive 
PRC), 2) the action potentials (spike widths) are narrow compared to their typical spiking period 
so the frequency is dominated by the membrane recovery time, and 3) the time scale for spike 
generation is very fast compared to the recovery time so an effective threshold for spiking can be 
defined. The first point was verified by observing that the measured PRC of the conductance-based 
models near the bifurcation to firing is positive. 

We model the approach to threshold with a simple passive decay to obtain dynamics governed 
by a single equation of the form: 

C^:^i-g^{V-Vr)-gS{i){V-Vs), (3) 
at 

where gm is an effective membrane recovery conductance, Vr is an effective membrane reversal 
potential and is a synaptic reversal potential. V{t) is reset to Vq whenever it reaches the threshold 
potential Vt- The passive decay to threshold is a very good approximation for some neuronal models 
such as the reduced Traub and Miles model given in the Appendix (Ermentrout and Kopell, 1998). 
On the other hand, we will show that even when the passive decay is not a good approximation 
to the slow dynamics of the neuron model, it can still adequately describe the frequency behavior, 
especially in the phasic regime where the synaptic current dominates. 

The synaptic current S{t) is generated from the spikes of pre-synaptic neurons. This must be 
emulated in the reduced model (|3|). Here we consider S{i) to be an arbitrary time dependent 
function. In Sec. ^, we analyze some biophysical synaptic models in detail and explicitly derive the 
time course of S{i) in response to a pre-synaptic spike. 

To simplify the analysis, we rescale Eq. (^) so that only dimensionless parameters remain. The 
voltage can be rescaled via v = {V — Vo)/{Vt — Vq), so that the reset potential is at f = and the 
threshold is at t; = 1. We define an intrinsic membrane decay time 

Tm = C/g.m (4) 

and rescale time by t = t/Tm- This scaling takes the membrane time scale to 1, leading to the 
dimensionless equation 

- = I-[l + ^S{t)]v-gS{t), (5) 

where 

^ I + Ir g g 

I = — , g = —, 7 = — , (6) 
It gr gm 
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with It = gmiVr — Vo), Ir = 9'm(K — Vo), and gT = /r/(Vb — Vs). The neuron is said to fire each 
time v{t) reaches 1, whereupon it immediately resets to v = 0. It is important to note that this 
does not restrict v{t) to nonnegative values away from t = 0. 
An even simpler model is 

^ = (7) 

which follows from Eq. (^) if 7 << 1. We refer to the reduced model given by Eq. (^) as model 
1 and the one given by Eq. (0) as model 2. The two models differ in that the synapse acts solely 
as a forcing function in model 2, but also affects the membrane decay rate in model 1. As we will 
show, the frequency behavior of the two models is qualitatively similar even if 7 is not small. We 
will focus most attention on model 2 in our analysis. The synchronization tendencies of model 2 for 
slow synapses has been studied thoroughly (van Vreeswijk et al, 1994; Hansel et al., 1995; Gerstner 
et al, 1996; Chow, 1998). 



3. Synaptic Model 

In conductance-based neuron models, the post-synaptic current is initiated by a pre-synaptic spike. 
In the reduced model, the spikes have been eliminated so this process must be modeled. We do so 
by deriving the synaptic time course for Eq. (^) explicitly for a stereotypical pre-synaptic spike. 
As we will show below, the spike response can be described by a time dependent recursive forcing 
function of the form 

S{t)^a{t-ti)S{t) + Sf{t-ti), (8) 

where S{t) is the synaptic gating variable to be used in Eq. (|5|), ti gives the spiking times of 
the pre-synaptic neuron, Sf{t) is the stereotypical post-synaptic response, and a{t) is a 'memory' 
function. (Note that all time and rate parameters have been rescaled by Tm in this section). We 
divide synapses into two types ~ saturating and nonsaturating. If the synaptic variable S{t) increases 
without bound as the rate of firing of the pre-synaptic cell approaches infinity then we call this 
a nonsaturating synapse. If however, S{t) saturates to a fixed value as the firing rate approaches 
infinity then we call this a saturating synapse. The behavior of the function a{t) determines the 
type of the synapse. 

We first derive the recursion relation explicitly for a simple nonsaturating synaptic model given 

by 

^ = -S/T + Y,S{t-ti), (9) 

1=0 

where ti denotes the times of the pre-synaptic spikes. Synapses of this form have been considered 
in previous models (Abbott and van Vreeswijk, 1993; Tsodyks et al, 1993; van Vreeswijk et al, 
1994; Hansel et al, 1995). Integrating Eq. (^ results in 

N 

S{t) = 5(0)e-*/" + e^^*"*'^/", (10) 
1=0 

where we begin the integration at t = 0. This sum is equivalent to the recursion relation 



S{t) ^ 5(t) + e-(*-*')/". 



(11) 



6 



Here a{t) = 1 and S{t) increases without bound as the frequency increases. 

Now consider the synaptic model in the Appendix. We will show that this is a saturating synapse. 
When the pre-synaptic neuron fires, -F(Vpre) in Eq. (^) rises from near zero to a value near unity, 
then returns to zero. The precise shape is determined by the temporal characteristics of the action 
potential and F{V). Our analysis is similar to that of Destexhe et al. (1994). We assume that the 
shape of F{V) can be approximated by a square pulse with a width of At, given by the time the 
membrane potential remains above zero during an action potential (i.e. width of the spike). We 
consider a train of square pulses arriving at times t;, with ti + At < t^+i VI (i.e. no overlapping 
spikes). The synaptic response is then found from 

^=aQ{t){l-S)-pS, (12) 

oo 

Q{t)= H{t-ti), (13) 

oo 

where H{t) is a square pulse of unit height with rise at t = and fall at t = At. 
Given /, we define 

Rit)^afQit')dt' = h^l-'^^^'^^'^'^^''' (14) 
^ ^ 4 ^ ' [aAt, ti + At<t<ti+i ^ > 

and then integrate (^) from ti io t < t^+i to yield 



S{t) = e-«(*)-'3* 



S{ti)e^'^+ /"*aQ(t')e^(*')+^*'dt' 



-a(t-t,)(5(t;)e-/3{t-t0) + (^_^^ (1 _ e-(a+/3)(t-t0) , ti<t<ti + At 

-^*(S(t/)e-/3(*-*')) + (^) (1 - e-("+^)^*) e-/3(*-*'-^*), ti + At<t< ti+i 



(15) 
(16) 



Consider S{t) at the time of the next spike, t = t^+i. If the spike does not occur, S{t) will continue 
to decay exponentially with rate P, i.e. 

S{t) = S{ti+i) exp(-/3(t - ti+i)), ti+i < t. (17) 

We can thus rewrite Eq. (|T^ ) as a recursion relationship with the spike at as 

S{t) ^ a{t - ti+i)S{t) + Sf{t - ti+i) (18) 

where 

r e-*, 0<t<At, 
= 1 e-^*. At < t, 

c.,^ 0<t<At, 
^n*) - I (_^) (1 _ e-("+/5)A*) e-/3(t-At)^ < t l^Uj 

The S{t) appearing in the right hand side of Eq. (|T8|) is understood to be the function that would 
occur in the absence of a new spike, Eq. (17). 

In contrast to Eq. (|ll]), there is a multiplicative factor a{t) which 'damps' away the past. For 
a{t) < 1, S{t) always saturates to a finite value. Sf{t) and a{t) have discontinuous first derivatives 
because a square pulse was used as input, but in general the functions will be smooth. 



and 
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Depending on the values of the three parameters a, (5 and Ai, the synaptic time course can 
take on many shapes. Here we are concerned with narrow spikes (small At) and slowly decaying 
inhibition (/3 << a). If the rise time is very fast, we can explicitly take the double limit a — > oo, 
At 0, a At = c, where c is a constant, to get 

5;(t) = (l-e-'=)e-^*. (21) 

and a = e~^. 

The constant c is the ratio of the spike width to the rise time of the synaptic current, and 
determines the contribution from the past to the current synaptic response. If c is large enough we 
can reasonably ignore the past and use the approximation a ~ and Sf{t) ~ exp(— /3t). Then the 
synaptic update function to be used in Eqs. ^ and (0) takes the form 

S{t) ^ e-(*-*')/", (22) 



where r = f/rm is the rescaled synaptic decay time. Equation (22) gives the recursion relation 
for the synapse where a(t) = and is valid in the limit where the rise time of the synaptic gating 
variable is very fast compared to the spike width. We note that a nonzero rise time is very important 
for synchronization (van Vreeswijk et aZ., 1994; Hansel et a/., 1995; Gerstner et at, 1996; Terman 
et al., 1997; Chow, 1998) but it is not as important for frequency determination. 



4. Period of a synchronized network 

In this section we will derive the form of the network period analytically for model 2 [Eq. (^] in 
terms of a transcendental relation for saturating and non-saturating synapses. We will also analyze 
the differences that would result if we considered model 1 [Eq. (^)]. In general the same qualitative 
behavior holds between the two reduced models. To calculate the period for the reduced models, 
we assume the neuron spikes repeatedly with a period T (the frequency is defined as / = T^^). We 
then integrate the membrane equation from the time of the last spike to the time of the next spike. 
The period is obtained self consistently by imposing the constraint v{T) = 1. 

Suppose the neuron last fired at t = 0, then for model 1 the membrane voltage obeys 

rT 

t;(r) = 1 = e-'^(^) / e^^''\l - gS{t')]dt' (23) 
Jo 

where (for model 1) 

^(t)= / [l + ^S{t')]dt'. (24) 

For saturating synapses, S{t) = exp(— t/r), which implies that 

/i(t) = t + 7(1 - e-*/^)r. (25) 

For model 2, this simplifies to /u(t) = t. 

Equation ( |2^ ) gives the membrane dynamics of the current spiking cycle in response to synaptic 
inputs from all the previous cycles. We can evaluate the combined synaptic input from all of the 
previous spikes if we assume periodic firing. For the neuron spiking in the past at times t = —IT, 
I = 0, 1, 2, . . ., we can re-express the recursion relation ^ as the sum 

oo 

S{t) =^a}Sf{t + lT). (26) 



For the saturating synaptic model, a = exp(— aAt) and Sf is given in Eq. This is a geometric 
series and can be summed to give 

S{t) = {l-a) (27) 
1 — ae ^ 1^ 

For model 2 (using ^{t) = t) we can integrate (^) explicitly using (^) to obtain 

Equation (^) is a transcendental relation that determines the spiking period and hence frequency 
of the synchronous network. If we can completely ignore the past (i.e. use synaptic model (^^) 
then a ~ and we obtain 

v{T) = 1 = J(l - e-^) - -^{e-^l^ - e-^). (29) 

r — 1 

For simplicity, we will conduct our analysis on this equation. 

For nonsaturating synapses, we can again calculate S{t) for periodic firing as we did for the 
previous case. Here, we find 

si^) = T^^- (30) 

For model 2 (using ^{t) = t) we can integrate (p3|) explicitly using (50) to obtain 



5. Frequency Regimes 

In this section we describe in detail the three asymptotic regimes discussed in the Introduction. We 
derive the parameter ranges in which each regime holds, and the dependence of the network period 
on the inhibitory decay time within each regime. 

5.1. Tonic Regime 

The Tonic regime occurs where 

T«l, T«r. (32) 



Note that unity in (32) corresponds to the intrinsic time scale of the membrane, which has been 
scaled away in our normalization, and r represents the ratio of the synaptic decay time to the 
membrane time scale. The relationship in (|32[) allows us to expand the exponentials in Eq. (p9D to 



linear order in T and T/r to obtain for saturating synapses 

T^{I-g)-\ (33) 
To derive the condition which must be satisfied for this regime we impose our original assumptions 



(32) on Eq. (pq). This leads to the condition 



{I - g) « min[T, 1] 



(34) 
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Condition (|34| ) can be considered to be the definition of the tonic regime. A strong apphed current 
(large /) and a weak synapse (small g) suffice to satisfy it. 

In the tonic regime, the synaptic decay time does not influence the period, as can be seen from 
Eq. (|33|). The synapse only affects the period through the average amount of inhibition the neuron 
receives over the course of one period. We see this immediately by replacing S{t) in Eq. (|2^) by its 
average {S{t)). After integrating, we obtain the relation 

l=.[/-5(5(t))](l-e-^), (35) 

where the time average of the synaptic input S{t) is defined as 

{S{t)) = ^£ Sit')dt'. (36) 
In the tonic regime to leading order this leads to the simple result of {S{t)) ~ 1 for saturating 



synapses and {S{t)) ~ r/T for nonsaturating synapses (using Eq. (31)). If we substitute {S{t)) into 



Eq. (^) and linearize in T and T/r (using (32)), we obtain the same period as given by Eq. (|3 



for saturating synapses. For nonsaturating synapses, we find the period obeys 

r~/-i(i + 5T). (37) 

Unlike the saturating case there is linear dependence on r. However, since g/I is small in this 
regime this dependence is small. 

The inclusion of the synaptic contribution to the decay rate in model 1 does not alter the results. 
In the tonic regime, T << r so we can expand Eq. (|25| ) to obtain 

fi{t) '^t + -ft = rt. (38) 

The inclusion of the synaptic contribution to the rate only changes the effective passive decay rate 
from 1 to r. Using this in Eq. (pS|) and expanding to linear order yields the the same period as 
given in Eq. (|3|). 



5.2. Phasic Regime 

The phasic regime occurs where 

1 « T, 1 « r. (39) 

This implies that T >> T/r which allows us to ignore with respect to e~^^'^ and 1. In this 
regime the neuron fires at a frequency which is slow but on the order of r~^. Applying condition 



(ID to Eq. (H) yields 

1 = 1- -^e-'^'\ (40) 



for saturating synapses. Solving for T gives 

9^ 



T = Tin 



(r-l)(/-l) 



(41) 



Equation ( |4l| ) shows that in this regime the period is approximately proportional to r and loga- 
rithmically dependent on the other parameters. We called this regime the phasic regime in White 
et al. (1997) because the synaptic variable S{t) is dominated by a phasic component with a period 



10 



proportional to r. To derive the condition on the parameters for the network to be in the phasic 
regime, we apply the condition T >> 1 to Eq. (p]). We then have the condition 



In 



g-r 



.(.-i)(/-i)J>>^" 

Thus, the phasic regime is obtained for g large compared to / — 1. We should note that I > 1 is 
the condition required for firing and thus must always be satisfied. 

We can compare condition ( ^2| ) to condition ( |3^ for the tonic regime. Both the phasic and the 
tonic regimes require fairly large r but it is the relationship between g and / that distinguishes the 
two regimes. The phasic regime is attained by a strong synapse and a weak applied current and 
the tonic regime is attained by the opposite. 

For nonsaturating synapses using Eq. ( |3l| ) the period takes the form 



r In 



9T + [T 

(r-l)(/-l) 



(43) 



The period is again proportional to r but the logarithmic factor is changed. 

The period is modified for model 1 but remains proportional to r. From Eq. ( [25| ) we can see 
that fi{t) is bounded by rt > fj.{t) > t, where r = 1 + 7 > 1. Hence the period is bounded between 
(|ll) and 

T = rln\- '-^ -1. (44) 

(rr — — r) \ 

Thus, T (XT but the proportionality constant is changed by a logarithmic factor. Here I > r must 
be satisfied for the neuron to fire. 



5.3. Fast Regime 

The phasic and tonic regimes are applicable for slow synapses, where the synaptic decay time is 
comparable to or slower than the period. Here, we investigate the fast regime where the synapse is 
very fast compared to the period; i.e. 

T«T. (45) 

The fast regime can occur for a range of r. 

We first consider r << 1 so that we can ignore e"^/"^ with respect to e~^. In this case Eq. ( [2^ ) 
for saturating synapses becomes 

1 -/(l-e""^) -5Te~^, (46) 
where we have also expanded to linear order in r. From this we obtain 

'gr + r 



r = In 



I -1 



(47) 



Note that the network period is a logarithmic function of the parameters. Thus, in this regime, 
the period is dominated by the membrane time constant, which is implicit in our scaling. For this 



regime to be valid we must impose condition (|45D on (|47|) . Thus we have the condition 

'gr + r 



In 



I -I 



» T. 



(48) 



This can be satisfied for / ~ 1 for arbitrary g or for g large enough and arbitrary / > 1. Note that 



if r is decreased from either the tonic or phasic regimes (see Eq. (|34D and Eq. (|42D), Eq. (^) is 
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still satisfied, implying that the network can enter the fast regime from either the phasic or tonic 
regimes. This occurs because the relationship between g and I distinguishes the latter two regimes, 
while the condition for the fast regime is primarily dependent on r. 

For r near to 1, Eq. ( p7[ ) is no longer valid since we cannot ignore e~'^l'^ with respect to e""^. 
By going back to (pg|) for r ~ 1 and T >> 1 we get the new period relation 

1 ~ /(I - e-^) - gTe~^ (49) 

which again can be satisfied (since T >> r ~ 1) for / ~ 1 and arbitrary g or for g large enough 
and arbitrary / > 1. 

The period is identical for nonsaturating synapses since the synapse is so fast it does not have 
a chance to reach saturation. The period also does not change much for model 1. This is because 
the synapse decays quickly (i.e. r << T). Thus over most of the period, /x(t) ~ f is a good 
approximation to Eq. (^). 



6. Relationship between Full and Reduced Models: Two Examples 

The three regimes analyzed above have been observed in numerical simulations of conductance- 
based models. The phasic and tonic regimes were reported in White et al. (1998) and the fast 
regime, for spiking and/or bursting neurons, has been observed previously (Friesen, 1994; Perkel 
and Mulloney, 1974; Skinner et al, 1994; van Vreeswijk et al, 1994; Wang and Rinzel, 1992). Here, 
we show that a stronger correspondence can be drawn between the reduced and the conductance- 
based models. We show that the period for the reduced model as a function of the parameters f, 
/, and g approximates the period obtained from simulations of the two conductance-based neurons 
given in the Appendix. Specifically we show that the period data T = T{I, g, f) obtained from the 
simulations of the conductance based models is well approximated by the period function for the 
reduced models obtained from the relation ([2^). 

The period relation (|28|) for the reduced model is in terms of dimensionless quantities T, I, g, 
and r (i.e. T = T{I,g,T)) . To compare to the numerical results we restore the dimensions using 
the scaling transformations given in (^, namely 

T= — , I=— , g= — . (50) 

Tm J^T QT 

Our task therefore is to find a single set of scaling parameters /r, It-, gx and Tm for which the 
reduced model period matches that of the conductance-based model. 

We first numerically generated (for each of the two models given in the appendix) three tables of 
period data. In the first table, we listed T as I varied for several fixed values of g and f, while the 
other two tables similarly listed T as g or f varied respectively. Each table was then considered a 
numerical function, T^^^f-^^I), Tf^j -^{g) or T^j -^(f), and we looked for a choice of scaling parameters 
( ^0| ) such that the reduced model period function T = T{I, g, r) implicit in (|2^) was closest to 
the numerical functions. We defined a measure of error as the sum of the absolute values of the 
maximal deviation between the reduced and conductance-based models (i.e. err = ^max \T — T\, 
where the maximum was taken over all period values within a table and the sum was taken over 
the three tables). The fit error between the conductance-based and the reduced model was then 
minimized using the function FMINS in the software package MATLAB. Critical to the success of 
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the fitting procedure was finding a good initial guess. This was obtained by adjusting the scaling 
parameters by hand and fitting the functions by eye. 

The memory function a = exp(Q!At) was fixed for each conductance-based model. At was defined 
as the spike width at 1/ = 0, while both models used a = 1. The results of the fits for both models 
are shown in Fig. Note that each neuron model has a single transformation set which is used 
in all three panels. The phasic and tonic regimes are evident in the figures. We did not explore 
the fast regime. The period derived from the reduced model (28) is seen to capture quite well 
the dependence of the period upon the parameters of the conductance-based models. Note that 
the White et al. model has a different spike shape from the reduced Traub and Miles model (see 
Fig. |l|). The latter has a significant negative overshoot while the former does not. Nevertheless, the 
quality of the fit is equally good in both cases. 

The fit tends to work best in the phasic regime. This is expected because the period behavior 
is dominated by the choice of synaptic model (which is the same in the analysis and simulations). 
In the tonic regime, the effect of the synapse is weak and the frequency characteristics are more 
like that of the uncoupled neuron. Thus differences in the intrinsic membrane dynamics will be 
the most prominent here. The period of the uncoupled reduced model obeys the classic result 
T = ln(//I — /*) (Tuckwell, 1989), and we showed in the tonic regime that the synapse affects the 
period only by changing the amount of tonic input the neuron receives. The period of the Type I 
neurons in the sense of Ermentrout (1996) behave as T = (/ — /*)~^/2 near the bifurcation point. 
Thus at least near the onset of firing, the frequency behavior as a function of / should not be the 
same. 

As the synapse becomes stronger, the influence of the synaptic dynamics begins to compete with 
the intrinsic membrane dynamics. Deep within the phasic regime the synaptic dynamics completely 
dominate and the period is primarily determined by the response of the neuron to inhibition. Our 
analysis showed that the intrinsic dynamics are not important in the phasic regime. However, the 
use of model 2 instead of model 1 for equation ( p8| ) will result in quantitative inaccuracies even 
well within in the phasic regime. 



7. Discussion and Conclusions 

In this paper, we show that a fully synchronized network can have its frequency determined by sev- 
eral parameters, including applied currents, synaptic strength and synaptic decay time. However, in 
regimes other than the phasic regime, in which parameters other than the synaptic decay dominate 
the network frequency, even mild heterogeneity can eliminate coherence (White et al., 1998). Indeed, 
in that work we showed that increasing the drive /, which increases the uncoupled frequency of 
each cell or the common frequency of a homogenous inhibitory network, can desynchronize a mildly 
heterogeneous network. This is true even if the percentage difference in natural frequencies over 
the population is held constant. Decreasing the drive sufficiently puts the network in a regime 
where inhibition suppresses spiking in less excitable cells. This also causes a loss of coherence albeit 
through a different mechanism. These two effects limit the possible frequencies at which the network 
can be synchronized, and the time scale of the synapse plays a crucial role in determining the width 
of this frequency band. 

With respect to the hippocampus, the analysis predicts the fast 200 Hz ripples observed in 
CAl (Ylinen et al., 1995), cannot be solely mediated by GABA^ synaptic inhibition, whose time 
scales are much slower than 5 ms. An additional mechanism such as gap junction mediated electrical 
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coupling may be responsible for the observed synchronous fast ripples. The 40 Hz gamma rhythm, 
which is on the order of the time scale of GABA^i, can be supported by inhibition alone. This was 
observed in experiments in hippocampal slices (Whittington et al, 1995; Traub et al, 1996a; Jefferys 
et al, 1996). Our results on the dependence of the period on the synaptic decay time constant 
clarifies the experimental observations. Our analysis predicts that for coherent oscillations, the 
network should be in the phasic regime. To compare with the analysis, we replotted the experimental 
data of Whittington et al. (1995) as period versus synaptic decay time in Fig. |^. We note that the 
period T is proportional to the decay time r, as predicted for the phasic regime. Furthermore, r/T 
has a value less than 1, consistent with the requirement for being in the phasic regime. 

For homogeneous networks, the time scale of the inhibition is not critical for the frequency 
of a synchronized network. A homogeneous network can stably synchronize at parameter regimes 
outside the phasic regime. However, for a mildly heterogeneous network, the synchronization mech- 
anism must not only draw together phases, but must also help create a common frequency. As 
shown in (Chow, 1998), inhibitory synapses in the phasic regime do provide such a mechanism 
for heterogeneous networks. We wish to point out that this frequency dependent synchronization 
is not necessary for all networks. For instance, consider a network coupled with (nonrectified) 
electrical coupling. In the fully synchronized state, the coupling currents disappear and hence do 
not contribute to the network frequency. With mild heterogeneity, the network still can maintain 
coherence, but there is no time constant to affect the network frequency. 

The connection of the frequency of a synchronized network with inhibition to the inhibitory decay 
time has previously been understood intuitively (Destexhe et al., 1993; Skinner et al., 1993; Kopell 
and LeMasson, 1994; Wang and Buzsaki, 1996). Here we present a simple example to clarify how 
all of the relevant time scales interact to produce the regimes in which that intuition is correct. 
These conclusions assume that the membrane potential between spikes is governed by a passive 
decay process that approaches a fixed threshold. The reduction from model 1 to model 2 simplifies 
the effect of the synapse. The accuracy of our analytical calculations depend on the validity of 
these approximations. Nevertheless, we expect the three regimes to still exist for a much larger 
class of conductance based models and models of synapses. Evidence for this conjecture was given 
in Section |6[ We note that Izhikevich and Hoppensteadt (1997) have shown that for weakly coupled 
systems an arbitrary Type I excitable neuron can be transformed into integrate-and-fire form by 
a piece- wise continuous change of variables. Our results show that this correspondence can hold 
beyond the weak coupling limit. 
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Appendix 

A. Neuron Dynamics 

For our physiologically based neuron we consider a single-compartment model with inhibitory 
synapses obeying first-order kinetics. The membrane potential obeys the current balance equation 

C^=i-lNa-lK-lL-Is, (51) 

at 

where C = l/^F/cm^, / is the applied current, I^a = QNaT^^hiVi — V]\ja) and Ik = gKn'^iVi — Vk) 
are the spike generating currents, II = QLiVi — Vl) is the leak current and Ig = gSj{t){Vi — Vg) 
is the synaptic current. 

The interneuron model in White et al. (1998) used parameters: gNa = 30 mS/cm^, gx = 
20 mS/cm2, gL = 0.1 mS/cm^, VNa = 45 mV, Vk = -80 mV, Vl = -60 mV, Vg = -75 mV. The 
activation variable m was assumed fast and substituted with its asymptotic value m = m^{v) = 
(1 + exp[— 0.08(?; + 26)])"^. The gating variables h and n obey 

dh _ hoo{v) - /i dn _ noo{v) - n 
dt Th{v) ' dt Tn{v) 

with h^{v) = (1 + exp[0.13(t; + 38)]y\ Th{v) = 0.6/(1 + cxph0.12(t; + 67)]), noc{v) = (1 + 
exp[-0.045(t; + W)])-\ and Tn{v) = 0.5 + 2.0/(1 + exp[0.045(t; - 50)]). 

The reduced Traub and Miles model (Ermentrout and Kopell, 1998; Traub et al, 1996b) used 
parameters: g^a = 100 mS/cm^, gx = 80 mS/cm^, g^ = 0.1 mS/cm^, VNa = 50 mV, Vk = 
-100 mV, Vl = -67 mV, Vg = -80 mV; m = moo(t') = amiv)/iam{v) + Pm{v)), where am{v) = 
0.32(54 + i;)/(l -exp(-(i; + 54)/4)) and /3^(i;) = 0.28(t; + 27)/(exp((i; + 27)/5) - 1); 

— = an{v){l - n) - Mv)n (53) 

with aniv) = 0m2{v + 52)/(l - exp(-(v + 52)/5)), Pniv) = 0.5exp(-(t; + 57)/40); h = hoo{v) = 
max[l - 1.25ra,0]. 

For both models, the synaptic gating variable Sj{t) is assumed to obey first order kinetics of the 
form 

^=aFiVj)il-Sj)-pSj, (54) 

where a = 1 and (3 are respectively the synaptic rise and decay rates and FiVi) = l/(H-exp[— V^]). 
Transmission delays are neglected. We are interested in the response of the system to changes in 
the applied current /, synaptic conductance g, and synaptic decay time f = /3~^. The ODEs were 
integrated using either a fourth-order Runge-Kutta method or the Gear method. 
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Figure 1. Example voltage traces of spikes for the A) White et al. (1998) interneuron and B) reduced Traub and 
Miles model of Ermentrout and Kopell (1998). 
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Figure 2. Shown are example slices of the numerically determined period for the conductance-based models (dashed 
lines) compared to the period of the reduced model obtained from Eq. ( |2^ ) (solid lines). Figures A), B) and C) 
show results for the White et al. (1998) model with memory coefficient a — 0.30, and scaling parameters Ir ~ 1.9155, 
It = 1.4337, = 12.0230, and gr = 0.0851. Figures D), E) and F) show resuhs for the reduced Traub and 
Miles neuron with a — 0.74, Ir — 1.3546, It = 1.6211, r„i — 16.1158, and gT = 0.1111. From top to bottom the 
slices have parameters: A) {g,f) = (2, 50), (1, 32.5), (1, 15), (0.05, 5) B) {i,g) = (1.64, 1), (5, 1), (5, 0.5), (5, 0.05), C) 
(/,f) = (1.74, 50), (1.74, 32.5), (1.74, 15), (5, 15), D) {g,f) = (2, 50), (1, 32), (1, 15), (0.5, 5) E) 
(1, g) = (1.085, 1), (4.825, 1.25), (4.825, 0.5), (4.825, 0.05) and F) (J, f ) = (2.02, 50), (2.02, 32.5), (2.02, 15), (4.825, 15). 
The phasic and tonic regimes are most clearly distinguished in B) and E) where the period is proportional to f in 
the upper curves and independent of r in the lower curves. 
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Figure 3. Plot of period vs. synaptic decay time for a hippocampal slice showing linear behavior indicative of the 
phasic regime. Responses were recorded from inhibitory neurons evoked by glutamate appUcation during a wash-in 
of 2 /iM pentobarbital (Whittington et al., 1995). Experimental data is provided courtesy of J. Jefferys. 



